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Abstract. We report on recent progress of the implementation of the similarity renormalization group 
(SRG) for three-body interactions in a one-dimensional, bosonic model system using the plane wave basis. 
We discuss our implementation of the flow equations and show results that confirm that results in the 
three-body sector remain unchanged by the transformation of the Hamiltonian. We also show how the 
SRG transformation decouples low- from high-momentum nodes in the three-body sector and therefore 
simplifies the numerical calculation of observables. 
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1 Introduction 

Renormalization group methods have become an impor- 
tant tool in modern physics [1 . In particular, for stud- 
ies of strongly-interacting many-body systems they fre- 
quently facilitate the correct calculation of observables. 
One of these methods is the similarity renormalization 
group (SRG), which has been extensively used in con- 
densed matter physics [2] but has recently become of im- 
portance also in nuclear physics [3l|4] . The SRG essentially 
constitutes a flow equation generating unitarily equiva- 
lent Hamiltonians, which after evolution possess certain 
features that usually make the calculation of observables 
easier. The form of the evolved potential will depend on 
the so-called generator that is an essential ingredient of 
the flow equation. Frequently used generators drive the 
two-body potential to the diagonal in momentum space. 
Off-diagonal elements are thus driven to zero and low mo- 
menta effectively decouple from large momenta. 

Such capabilities are very important in ab initio nu- 
clear structure physics for which high-performance com- 
puting has become of increasing importance and the avail- 
able computational resources set a hard limit on the num- 
ber of observables that can be calculated [5 . Recently, 
the SRG flow equations for three-body interactions were 
implemented for the truncated harmonic oscillator basis 
and first nuclear structure calculations with consistently 
evolved two- and three-body interactions have emerged [6j 
[TllH] • The SRG therefore seems to provide a way to extend 
the limits of possible computations. In these calculations 
the interactions were evolved in the harmonic oscillator 
basis. However, a calculation of the evolved three-body 
interaction in the plane wave basis would not only fa- 
cilitate the projection of the potential on any basis but 
also the calculation of infinite matter observables such 



as the energy per particle of nuclear matter. In Refs. [9l 
[To] it was shown that an evolved two-body interaction in 
combination with a phenomenological three-body interac- 
tion seems to lead to a converging many-body perturba- 
tion theory series for infinite matter. However, a missing 
link in this calculation remains the inclusion of a consis- 
tently evolved three-body interaction. Instead, parameters 
of the leading chiral effective field theory three-nucleon 
force were refitted and used alongside the evolved two- 
body potential. Since these positive results for light nuclei 
and nuclear matter validate also an effort to construct an 
ab initio density functional based on many-body perturba- 
tion theory with evolved interactions [11] , the SRG seems 
therefore to have impact on calculations across the whole 
chart of nuclides. 

In this work we consider the evolution of a one-dimen- 
sional Hamiltonian with two- and three-body interaction 
terms. In Sec. |2]we will introduce the flow equations that 
determine the evolution of the Hamiltonian. In the follow- 
ing two sections we will then discuss how the flow equa- 
tions are implemented for two- and three-body interac- 
tions in bosonic systems. The presentation of numerical 
results will focus on illustrating key properties of evolved 
interactions, i.e. the conservation of observables and (for 
our choice of generator) the effective decoupling of low 
and high momenta. In the last section we summarize our 
results and discuss necessary steps towards an extension 
of this work to three dimensions and to nuclear systems. 



2 The Similarity Renormalization Group 

Let us denote with Hs the transformed (and initially un- 
known) Hamiltonian where s denotes the so-called flow 
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parameter and provides a measure of how much the Hamil- 
tonian has been transformed. Then, there exists a unitary 
transformation U. such that 



H, = UsHUl 



(1) 



where H = Hg=o is the known, original Hamiltonian. We 
can now calculate the derivative of the above equation 
with respect to s 



dH, dU, 



ds 



ds 



HUl + UsH 



dUl 
ds 



(2) 



Since the transformation Us is unitary we have UgUl = 1 
and therefore 



ds ^« " ^' ds 
Using this in Eq. ^ gives 
dH. 



ds 



rjsHs - HgTis = [r]s,Hs\. 



(3) 



(4) 



We can then specify the unitary transformation by speci- 
fying 7^5, which is subject to the condition 



V. 



t 



-Vs 



(5) 



which follows from Eq. ([3|. 

A convenient choice of 77^ is 77^ = \Gs-,Hs\ where Gs 
is a Hermitian operator. It obeys Eq. (pi) since Hg is also 
Hermitian: 

[Gs,Hs]^ = -[Gs,Hs]. (6) 

We will store all dependence on the flow parameter s in 
the potential term of the Hamiltonian, thus writing Hs = 
Trei + T4, where T^ei is the relative kinetic energy operator. 
Expanding the commutators then gives the equation 



dH, dK 



ds 



ds 



= GsHgHs + HgHsGs — 2HsGsHs. (7) 



There is significant freedom in the choice of the generator 
Gs- In this work we have chosen Gs = T^ei- This generator 
is used since it is known to drive the Hamiltonian to the 
diagonal, which is usually desirable. The reason for this 
feature is that T^ei is, in itself, diagonal in momentum 
space. 



3 Momentum Space Equations 

We will define the interactions with respect to the usual 
three-body Jacobi momenta 



q=^fk3-i(ki+k2)) . (8) 



We will be working in a partial-wave projected basis. This 
needs to be clarified since in one dimension plane waves 



can only propagate in two directions; backwards and for- 
ward. In one dimension, there exist therefore only two 
partial waves / = and / = 1, which correspond also to 
the parity of the state. We will be interested in bosons in 
this work and therefore assume that the two-body system 
is symmetric under exchange of the particles. Thus, we 
will only work with / = basis states. In the two-body 
sector, a complete set of states will therefore be written 
as 



J^dp \pl = 0){pl = 0| = J dp \pO){pO\ 



(9) 



In the three-body sector we will concentrate on states with 
total angular momentum L = 0, which implies that the 
relative angular momentum associated with the q variable 
is A = 0. The complete set of states in the three-body 
sector is therefore 



1 



dp / dq \pq (00)0) {pq{00)0\ 



(10) 



where the parentheses denote that two angular momenta 
have been coupled to total angular momentum L. From 
now on we will drop all angular momentum information in 
the bras and kets since we have only one angular momen- 
tum state in the two- and three-body system, respectively. 
Another quantity useful to define is the hypermomen- 
tum (^, defined by 

e=p' + lq'- (11) 

Since this quantity is proportional to the total, relative 
kinetic energy of the three-body system it also defines a 
plane in which the three-body interaction will become di- 
agonal through evolution. 



3.1 Interactions 

Two- Body Potentials: We have used different two-body 
potentials in this work to analyze the features of the SRG 
evolution. In particular, we have used a separable poten- 
tial that has the advantage that the binding energy and 
the two-body t-matrix can be calculated analytically 

Vsep{p,p') = 5exp( VM') exp(-p'V^') • (12) 

As the regulator A is increased towards infinity the poten- 
tial takes the form of a delta- function in coordinate space. 
In this case the binding energy of the A/'-boson state is 
known analytically and provides an excellent test for our 
few-body code. For additional benchmarking of numerical 
results we have also employed a potential that was previ- 
ously used in Ref. [12 



V{p,p') = X] ^expf 



ip-p') 



V? 



(13) 



with parameters given in Table [T] It is important to note 
that we employ the partial-wave projected versions of the 
above potentials. The parameters in Table [l] show that 
we will use a purely attractive interaction (Vg) and an 
attractive interaction with short-range repulsion (Va). 
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Table 1. Parameters for the two-body potential given by 
Eq. p). 





Vi 


V2 


cri 


cr2 


v^ 

Vp 


12 



-12 
-2 


0.2 



0.8 
0.8 



The Three- Body Potential: The SRG evolution will gen- 
erally induce many-body forces, however, it is expected 
that the three-body force will dominate over higher many- 
body forces as long as the flow parameter s is not too large. 
To complete our analysis, we also added a three-body force 
to the unevolved Hamiltonian to mimic general features 
present in nuclear physics. The three-body potential that 
we have used is of the same simple form as the one used 
in Ref. [12] 



where ce is the strength of the interaction and 



(14) 



fA{p,q) =exp 






(15) 



where we use n = 4 and A = 2 throughout this work. 



3.2 SRG Equations 

Two- Body SRG: With the kinetic energy operator as gen- 
erator for the SRG equation, Eq. pj becomes 



d^ 
ds 



2TVsT + V,V,T + TV,V,-VsTT-TTVs-2VsTV,. 

^ (16) 
Using the complete set of states defined in Eq. Q we can 
write out the evolution equations in momentum space 

±{p\V„\p') = - {p^ - p'^f {p\V.s\p') 



+ {p'+p'')J dq{p\VM){q\Vs\p') 

-2J dqq^{p\V,\q){q\Vs\p') ■ (17) 

The first term on the right hand side ensures that the 
potential is driven to the diagonal. 



Three-Body SRG: The equations for the evolution of the 
three-body potential are significantly more complicated. 
We write the Hamiltonian in the three-body sector as 



Hs=T^V, 



(1) 



V, 



(2) 



■Vr 



(3) 



^3, 



(18) 



rii) 



where V2 denotes the two-body potential in the three 
different channels. We have also dropped the subscript s 
from the potential and will keep doing this from now on 
to simplify notation. Since we are working with identical 



bosons we will assume that the induced three-body inter- 
action is symmetric under the exchange of two particles [^ 
The flow equation in the three-body sector is then written 



as 



liv^^' + vP 



V^^^+Vs) = [[T,H,],H,]. (19) 



The expression above contains the combined evolution of 
the two- and three-body interaction. It is of general inter- 
est to separate the two-body from the three-body evolu- 
tion but in this case it removes also spectator (^-functions 
that arise from disconnected diagrams and complicate the 
numerical computation of the evolved potential. Follow- 
ing [3 , we circumvent this issue by subtracting from the 
above expression the evolution of the two-body potentials 
V2 ' This isolates the derivative of the three-body poten- 
tial and removes aforementioned (5-functions 



ds 



[[T^Hs],Hs, 






ds 



(20) 



Expanding the commutators and rewriting the two-body 
differential equations as in previous section gives us 



ds 
where we have defined 
3 



O2 + 02^ + O^ , 



(21) 



02= E (^ - ^^0) (tV^^V^'^ + V^W^'^T - 2V^^TV^'^) 



3 



023 = ^ TV^V^'^ + V^V^^T - 2V3TV, 



(i) 



O3 = 2TV3T - 2V3TVS + 

TV3V3 + V3V3T-V3TT-TTV3. (22) 

(2) (3) 

We can express the two-body potentials V2 and V2 

through the potential V2 after application of permuta- 
tion operators 



V, 



Pl2,P22y2 P12P22, 

-(3) _ D D„ T/(i); 



V^""^ = Pl2P23V2''^Pl3P23 • (23) 

At this point it is useful to define the operator P given by 
P = P12P23 + P13P23 , (24) 

where Pij denotes the permutation operator that exchanges 
particles i and j. It can be shown that the overlap ma- 
trix elements for P12P23 and P13P23 in a partial-wave pro- 
jected basis are identical. We can therefore write 



V, 



(2) 



1 



-(l)i 



V, 



(3) 



^PV2'''P , 
]PV2'''P , 



(25) 



^ In the fermionic case we will write V3 = Vo + Vo + Vo 
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which simphfies the above equations significantly. 

We have used two different representations of the ma- 
trix element of the operator P 



s = 



s = 0.001 



s = 0.01 



{pq\P\p'q') = "^Sip- 7r((7, q' , x))5{p' - 7r{q\ q, x)), 

x=±l 

(pq\P\p'q') = E ^(P - ^(^'' «'' ^))^(« - ^(^'' «'' ^))' 

(26) 



X^±l 



where 



7r((7, q, x) = J -q^ + q''^ + xqq' , 



Tf{p, q, x) = \ -p2 



16' 



-pqx, 



X(p,^,x) = \Ip'^ + -q^ -pqx. 



(27) 



The obvious consequence of the implementation of the 
operator P is off- grid momenta in the object it is applied 
on. This problem can be solved by splining these objects, 
e.g. a function containing the shifted momentum 7r(g', q' ^ x) 
will be written as 



N 



f{7r{q,q\x)) = "^Siiniq.q' ,x))f{qi) 



(28) 




■M 



10 20 10 20 10 20 

P P P 

s = s = 0.001 s = 0.01 




20 



20 



s = 



s = 0.001 



s = 0.01 




20 



20 



We have used the global splines defined in Ref. [16 but 
also the cubic splines given in Ref. [17]. While the cubic 
splines provide a speedup in the calculation they also de- 
crease the accuracy of results slightly. The results in this 
work were therefore all obtained with global splines. 



4 Observables 

4.1 Two-Body Observables 

Phaseshifts and binding energies: We have calculated 
scattering and bound-state properties in the two-body 
sector. Scattering properties are obtained by solving the 
Lippmann-Schwinger equation. In operator form it is given 
by 

t = V^VGo{E)t , (29) 

where Gq{E) denotes the free Green's function 

1 



Fig. 1. SRG evolution of the potentials Va (upper panel), Vp 
(middle panel) and 14ep (lower panel) as defined in Eqs. (13) 
and (12) and Table fl] The parameters in the separable poten- 
tial are given by ^ = — 1 and A — 10. 

Table 2. Binding energy B2 for two particles interacting 
through the evolved potential \4ep for different s. The origi- 
nal potential is given by Eq. (12) with g — —1 and A — 10. 



s 


B2 





5.12880487 


0.005 


5.12880468 


0.01 


5.12880466 


0.02 


5.12880462 


0.04 


5.12880457 


0.06 


5.12880454 



Go{E) 



E — p^ jm + is 



(30) 



The phaseshifts ^i are then obtained from the on-shell t- 
matrix using the relation 



4.2 Three-Body Observables 

In the three-body sector we will focus on binding energies. 
An effective way to calculate three-body binding energies 
in momentum space is provided by the Faddeev equation. 

IV) = Got2P|V> + Got2Got3(l + Pm. (32) 



tiiP,p) 



2pe^^' sin 5i 



(31) 



where ^2 and ts are the transition operators ( 29 ) obtained 



where we will concern ourselves only with the / = phase- 
shifts. 



from the two- and three-body potential terms. 
We can then obtain the total wavefunction by 

\^) = {i+Gots){i+pm. 



(33) 
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Fig. 2. Two-body scattering phaseshifts calculated from the 
evolved and unevolved potential V^- The different panels cor- 
respond to different values of the evolution parameters s with 
s = (upper panel), s — 0.001 (middle panel) and s — 0.01 
(lower panel) . The different curves denote results obtained with 
different values of the Lippmann-Schwinger momentum space 
cutoff vlcut- The solid line denotes the uncut result. The dashed, 
dot-dashed, dotted lines give the result for ylcut = 7, 5, 3, re- 
spectively. 



Table 3. Two- and three-body binding energies B2 and B^ 
for the starting two-body potentials Va and Vp and varying 
strength of the s = three-body interaction. 



V2 


CE 


B2 


B3 


Vo. 


-0.10 


-0.920 


-3.226 


K 


-0.05 


-0.920 


-2.885 


K 


0.00 


-0.920 


-2.567 


K 


0.05 


-0.920 


-2.279 


Vc 


0.10 


-0.920 


-2.027 


Vp 


-0.10 


-0.474 


-2.570 


V, 


-0.05 


-0.474 


-2.132 


y& 


0.00 


-0.474 


-1.708 


y& 


0.05 


-0.474 


-1.307 


Vp 


0.10 


-0.474 


-0.952 



programming language. We did not encounter any stiff- 
ness when solving the differential equations. 

We have chosen to work in units where h^ /m = 1 
and all results are given according to this convention. The 
results for two-body binding energies obtained from un- 
evolved potentials are shown in Table [3] and the first row 
of Table [2] Three-body binding energies are presented in 
Table [3J We note that our results for V5 + V3 differ slightly 
from tne values given in Table II of Ref. [12]. However, our 
results have been reproduced recently \JE and are there- 
fore assumed to be correct. 




Fig. 3. Two-body binding energy of the evolved poten- 
tial Vc, as defined in Eq. ( 13 ) and Table fl] The y-axis de- 
notes the momentum-space cutoff ylcut in the momentum-space 
Schrodinger equation. The different lines correspond to differ- 
ent values of the evolution parameter s with s = (solid line), 
s = 0.001 (dashed line), s = 0.01 (dot-dashed) and s = 0.1 
(dotted line). 



5 Results 

We have implemented the evolution of the two- and three- 
body potentials in Python and Matlab. Equations (17) 
and (21) were discretized and written as a number of ma- 
They were then solved using one of 



trix multiplications. 

the standard ode solvers available in the corresponding 



Two- Body Evolution/Decoupling: We have evolved all 
previously defined two-body interactions. The diagonal- 
ization as a result of the evolution in the plane that is 
spanned by incoming and outgoing relative momenta is 
common to all starting interactions. In Fig. [l] we show 
the SRG transformation of interactions V^^ Vf3 and l^ep, 
respectively. With increasing flow-parameter s the poten- 
tial becomes increasingly more diagonal. The area in the 
low-momentum region of these figures that does not get 
diagonal with the evolution indicates the presence of a 
low-momentum two-body bound state. This is a typical 
feature of the T^ei generator and is not be present for all 
other possible generators such as the Wegner generator 
used in Ref. [13 . 

The decoupling of large and small momenta in the 
potential through the SRG evolution is one of the most 
important features that result from employing the diago- 
nal Trei generator. Low-energy observables can therefore 
be calculated correctly with a decreased momentum-space 
cutoff after sufficient evolution. We illustrate this impor- 
tant feature of the SRG in Figs. |2] and [3] where we plot the 
phaseshifts and two-body binding energy (obtained with 
V^) as a function of a sharp cutoff ylcut in the momentum- 
space Schrodinger equation for different values of s. With 
increasing flow parameter 5, observables become less sensi- 
tive to such a truncation indicating thereby the decoupling 
of small from large momenta in the potential. 



Three- Body Evolution/Decoupling: We have evolved all 
potentials and obtained not only the corresponding two- 
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body potential but also the induced three-body interac- 
tion term. Our generator diagonalizes the Hamiltonian in 
terms of the total kinetic energy. In Figs. |4] and [5) we have 
therefore chosen to plot the three-body interaction in the 
plane of incoming and outgoing hypermomentum (" and 
C,' for three different hyperangles = arctan [\/?>p / {2q)) . 
The induced three-body potential obtained in the absence 
of an initial three-body force is shown in Fig.|4] Note that, 
in this case, the three-body potential is identical to zero for 
5 = 0. It can be seen clearly how the evolution induces a 
three-body potential and that the strength of off-diagonal 
elements in hypermomentum space depends strongly on 
the flow parameter s. The three-body potential is weak for 
small s parameter but also couples strongly off-diagonal 
elements in the hypermomentum plane. Then, when s in- 
creases the three-body potential becomes stronger but also 
more diagonal in the hypermomentum plane. In Fig. [5] 
we show the same evolution of the three-body interac- 
tion where the Hamiltonian contains the initial three-body 
force as defined in Eq. (14) with ce = —0.05. 



Furthermore, we have tested decoupling in the three- 
body sector numerically. This was achieved by multiplying 
the two- and three-body potentials with regulators that 
are functions of the two- and three-body kinetic energy, 
respectively. In particular, we modify the evolved potential 

V,{p,q,p',q') ^ exp(-(C' + C^)/Al,,)Vs{p,q,p',q') . 

(34) 
In Fig. [6] we show the three-body binding energy as a 
function of ylcut foi" different values of the flow parame- 
ter s. The longer evolution, measured by s, the smaller 
is the minimal truncation cutoff ylcut for which the three- 
body binding energy remains unchanged. The decoupling 
of small and large hypermomenta due to the SRG is there- 
fore clearly visible. 



6 Summary and Outlook 

In this work we have presented results for the SRG evolu- 
tion of three-body interactions for one-dimensional, bosonic 
model systems performed in a plane- wave basis. We showed 
that the evolution in the plane-wave basis is easily im- 
plemented and that observables in the three-body sector 
remain unchanged when evaluated with the evolved po- 
tentials. We showed explicitly that the SRG induces a 
three-body potential that becomes increasingly diagonal 
in the hypermomentum plane when evolved to larger flow 
parameter s. We tested decoupling explicitly by calculat- 
ing observables with truncated two- and three-body inter- 
actions. Decoupling works well as illustrated for binding 
energies and phaseshifts in the two-body sector and for 
binding energies in the three-body sector. 

This work presents the first step towards a consistent 
SRG evolution of two- and three-nucleon interactions in 
the plane- wave basis. The extension to three dimensions 
and to nuclear systems will require small modifications 
of the spatial part of the equations plus the addition of 
angular momentum, spin and isospin recoupling to the 
problem. The latter is a usual part of standard few-body 
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Fig. 4. The induced three-body potential as a function of the 
inital- and final-state hypermomentum (" and C for the hy- 
perangles = 7r/12 (upper panel), = 7r/4 (middle panel), 
= 7r/2 (lower panel). The starting two-body potential is Va, 
as defined in Eq. ([l3} and Table [l] 



equations and we do therefore expect these changes to be 
straightforward. The immediate benefit of this implemen- 
tation is the increased consistency in a calculation of nu- 
clear matter observables. A plane- wave formulation will 
also allow to test more reliably the importance of four- 
body forces in a description of nuclear systems. It was 
claimed recently that these become of increasing impor- 
tance (after SRG evolution) in certain nuclei [8j. However, 
the results in Ref. [8 were obtained with a SRG evolu- 
tion in the harmonic-oscillator basis that could in prin- 
ciple be associated with truncation errors. Furthermore, 
a formulation in a plane-wave basis opens up the possi- 
bility of identifying relevant scattering observables in the 
four-body sector where an emerging four-body interaction 
should also be visible. 

Another possible avenue for the SRG in three dimen- 
sions is the Efimov effect in three-body systems of iden- 
tical bosons with large scattering length. It is an open 
question how the renormalization group limit cycle, that 
was found to occur in this problem [18], manifests itself 
as the SRG parameter is varied and the induced three- 
body interaction changes. Decoupling in the three-body 
sector might also provide another path to extract univer- 
sal properties of the four-body sector. The SRG might be 
used to remove three-body bound states and to calculate 
the universal properties of highly-excited four-body states 
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Fig. 5. The evolved three-body potential as a function of the 
inital- and final-state hypermomentum C, and C,' for the hy- 
perangles = 7r/12 (upper panel), = 7r/4 (middle panel), 
— 7r/2 (lower panel). The starting two-body potential is Va 
as defined in Eq. ( 13 ) and Table [l] the starting three-body 
interaction is given in Eq. (14), where ce — —0.05. 



without the need of a scattering calculation as done for ex- 
ample in Ref. [19]. 
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